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Abstract: We detected the H 2 6.3 /im bands in more than 30 normal red giants stars 
from K5III to M8III as well as in some early M supergiants on the SWS spectra retrieved 
from the ISO Data Archive. This result, however, shows serious inconsistency with the 
present model photospheres which predict H 2 only in the latest M (super)giant stars. 



Also H 2 was once discovered in the early M (super)giant stars nearly 40 years ago with 
psj \ the balloon-borne telescope named Stratoscope II. This discovery was so unexpected at 
that time that it was not understood correctly and overlooked for a long time. Now, we 
reflect on our ignorance of this important discovery during the 40 years and should consider 
more seriously the meaning of the rediscovery of water in so many red (super)giant stars 



P i with ISO. 
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1. INTRODUCTION 

At this monumental epoch to celebrate the opening of the ISO active archive phase, it 
may be instructive to recall a short history of stellar spectroscopy in space. At the infancy 
of the infrared astronomy in the 1960 's, an ambitious attempt to observe stellar spectra 
with a balloon-born telescope was undertaken and this mission named Stratoscope II, 
launched on March 1963, successfully observed infrared spectra (0.8 - 3.1 fim) of several 
red giant and supergiant stars (and invaluable spectra of Jupiter). The results showed 
beautiful spectra of water in Mira variables o Cet and R Leo (Woolf et al. 1964). This 
result was well in accord with the theoretical prediction (Russell 1934) and thus was well 
appreciated at that time. 
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However, the Stratoscope observers reported a more surprising result that water was 
detected in the earlier M giants [i Gem (M3III) and p Per (M4II) as well as in the early 
M supergiants a Ori (M2Iab) (Woolf et al. 1964) and n Cep (M2Ia) (Danielson et al. 
1965). However, this discovery was not in accord with the understanding of cool stellar 
atmosphere at that time and it was finally reinterpreted that the absorption bands at 
0.9, 1.1, 1.4, & 1.9 /xm identified with H 2 by the Stratoscope observers should instead 
be due to the CN Red System which also has band heads at about the same positions 
(Wing & Spinrad 1970). This proposition was more easily accepted by the astronomical 
community since then, since CN had been observed in a wide range of oxygen-rich stars 
from the Sun to red supergiants, not to speak of carbon stars. 

Meanwhile, the discovery of water, at least in the coolest Miras, confirmed the impor- 
tance of water as a source of opacity in cool stars, and actual computation of the non-grey 
model photospheres revealed that this is true in red giant stars with T eS cooler than about 
3200 K (Tsuji 1978), which roughly corresponds to M6III (Ridgway et al. 1980). This 
result was well consistent with the detection of the strong H2O bands in the Miras but 
was contradicting with the identification of H2O in the M2 - M4 (super)giants by the 
Stratoscope II observers. Thus, this result lent further support to the Wing-Spinrad 
proposition. A blind spot in this apparently reasonable conclusion, however, was that we 
were not aware at that time that the model photosphere, often referred to inadvertently 
as model atmosphere, is simply a model of the photosphere (0 <^ t < 00) and not of the 
atmosphere (—00 < r < 00), which may still involve unknown problems. 

Since then, important infrared missions such as IRAS and COBE have been success- 
fully undertaken, but few observations of the near infrared stellar spectra were done in 
space, except for continued efforts with KAO which provided fine low resolution stellar 
spectra (e.g. Strecker et al. 1979). Finally, ISO launched on November 1995 (Kessler et 
al. 1996) provided the means by which to observe astronomical spectra at higher resolu- 
tions for a wider spectral coverage at last. One of the unexpected results in the initial ISO 
observations with SWS (de Graauw et al. 1996) was a detection of water in the M2 giant 
f3 Peg (Tsuji et al. 1997) and in several early M supergiants in h + x P er clusters (Tsuji 
et al. 1998). It took sometime before we recognized that the Stratoscope II observers 
correctly identified water in the M2 - M4 (super)giant stars and that the Wing-Spinrad 
proposition was not correct (Tsuji 2000a). Further, IRTS launched by ISAS on March 
1995 detected H 2 bands at 1.9 fim in several M (super)giants earlier than M6 (Matsuura 
et al. 1999). Also, ground-based mid-infrared spectroscopy revealed H 2 pure- rotation 
lines in a Ori and a Sco (Jennings & Sada 1998). 

All these results that H 2 exists in red (super)giants earlier than about M6, however, 
could not be understood with the present model photospheres. We now encounter a 
serious problem: are we confronting with a fall of the classical model stellar photosphere 
or else with a rise of a novel picture of the stellar atmosphere (including the photosphere 
as a part)? A rather intriguing story of water in red (super)giant stars still remains 
open, and we hope to utilize the extensive ISO Data Archive to extend and finalize this 
fascinating story. 



2. A RED GIANT SAMPLE 



Our initial detection of water was done with the H 2 2.7 fim bands which, however, are 
contaminated by OH, CO, and other molecular bands. Then, we analyzed the H 2 6.3 fim 
bands which are little disturbed by other molecular bands. The expected spectra of the 
H2O V2 bands computed with the use of HITEMP (Rothman 1997) based on a single 
absorption slab model are shown in Fig.l. 

We first examined a dozen of high resolution spectra of red giants in the ISO Data 
Archive. The H 2 6.3 fim bands are detected in the K5 giant Aldebaran (but not in other 
two K giants) as well as in all the M giants between MOIII and M3.5III (Tsuji 2001), 
and also in the later M giants (M6-7III), as shown in Fig. 2. We found T ex « 1500K (and 
log iV co i noted on Fig. 2) by referring to the spectra such as shown in Fig.l. It is remarkable 
that SWS detected such faint water bands in the K giant star a Tau. This result is quite 
unexpected but confirmed recently in another K giant Arcturus (a Boo) with the high 
resolution ground-based spectroscopy in the 12 ^m region (Ryde et al. 2002). 

Next, we extend our survey to a larger sample of the low resolution SWS spectra. 
We found dozens of spectra in this category from the ISO Data Archive. At the lower 
resolution, some details of the band structure seen at the higher resolution are smeared 
out (Fig.l). Nevertheless we can detect the dip at 6.63 fim due to the H 2 is 2 bands in 
25 M giants from MOIII to M8III. Some examples are shown in Fig.3 and we estimated 
iV co i values for these M giants, again with the single slab model of Fig.l. 

The resulting values of N co \ from the high and low resolution samples are plotted 
against spectral types in Fig.4. For comparison, the predicted values of iV co i from the 
spherically extended non-grey model photospheres are shown by the dashed line. It is 
clear that the observed iV co i values cannot be explained at all by the predicted ones. 
Thus, H 2 detected in M giants should be non-photospheric in origin, but where does it 
come from? 



3. A RED SUPERGIANT SAMPLE 

In red supergiant stars, the H 2 6.3 ^m bands are found as absorption in a Ori 
(M2Iab) and a Sco (M2Ib), as in K - M giants. Further, water appears in absorption at 
A < 5 fim but in emission at A ;> 5 fim throughout in fi Cep (M2Ia) (Tsuji 2000b). This 
detection of water in distinct emission should be an important clue to the origin of H 2 0, 
since such emission should most probably originate in the outer atmosphere and not, for 
example, in the "starspots" . 

To account for the emission, we upgrade our single slab model to a spherically extended 
molecular sphere (MOLsphere for simplicity). Since we already know that T cx « 1500 K 
and N co \ 3 x 10 20 cur 2 from the Stratoscope II data (Tsuji 2000a), an additional free 
parameter is the inner radius of the MOLsphere r^. For simplicity, we consider only H 2 
whose absorption cross-section is as large as 10 -18 cm 2 and thus the H 2 gas is optically 
thick. Then we solve the transfer equation with the photospheric radiation, resulting in 



Fp, as a boundary condition, which shows absorption bands due to CO, CN, OH, SiO 
etc (Fig. 5). We found that the resulting emergent flux from the MOLsphere F P+M with 
r\ fa (R* is the stellar radius) accounts for the prominent emission lines due to H 2 
at A ^ 5 /xm as well as the absorption bands at A < 5 (Tsuji 2000b). 

Further, we found that the huge infrared excess can be explained by an optically thin 
dust envelope with rf ~ 0.1 and rf ~ 13. 5R*. Then the dust emission is simply 
added to F P+M - The entire spectra of [i Cep observed with the ISO SWS (corrected 
for Ay = 1.5 mag.) can reasonably be explained by our final spectrum resulting from 
Photospher + MOLsphere + Dust Envelope, Fp+m+d- Some details of the water emission 
in the 6 jum and 40 jum regions can be well reproduced by our model as shown in the 
inserted boxes of the left and right, respectively, in Fig. 5. So far, the huge IR excess of 
H Cep was thought to be due to dust alone, but now it is clear that it includes water 
emission originating in the MOLsphere. 



4. DISCUSSION 

Although we have assumed MOLsphere for fi Cep, this is not a theoretical model of the 
usual meaning but simply a kind of working hypothesis or an empirical model at best. 
We know nothing about the origin of the MOLsphere and it is a major challenge how 
to resolve this issue. Such a difficulty, however, may be shared with the origin of the 
chromosphere as well as of the mass- loss outflow (at least in non-pulsating stars). With 
this reservation, our empirical model for n Cep is reasonably successful, and a problem 
is if such a model can be extended to other cases. Although we see no emission in the 
6 fxm region in our large sample of red giant stars, it is interesting to notice that the H 2 
column densities tend to lever-off at about M5 and then to be smaller than the predicted 
photospheric values (Fig. 4). Moreover, water bands almost disappear in M7 giant EP 
Aqr. These results suggest that the photospheric H 2 bands may be filled in by the 
emission due to H 2 itself in the late M giants. Also, in the late M giants, H 2 appears 
as emission in the 40 ^m region (Tsuji et al. 1999), C0 2 shows prominent emission in 
the 15 ^m region (Jasttanont et al. 1998), and S0 2 shows emission as well as absorption 
in the 7^m region (Yamamura et al. 1999). Based on these observations, it should be 
reasonable to assume the presence of MOLsphere in late M giants. 

However, a possibility that there are some serious flaws in the present model photo- 
spheres cannot be excluded, especially for K (see Ryde et al. 2002) and early M giants for 
which there is no direct evidence for MOLsphere. This case offers a more serious problem, 
since our present stellar spectroscopy (e.g. abundance analysis) loses its basis. Also, if 
this is the case, we must give up to have a unified picture for cool luminous stars including 
K and early M giants. We are yet tempted to have a unified picture, ISO's innovation of 
which may be illustrated as in Fig. 6. The presence of the high excited molecular gas has 
been known for late (super)giants from the H 2 and SiO masers for a long time, but ISO 
revealed that the presence of such a warm molecular gas, including not only H 2 and SiO 
but also C0 2 and S0 2 , may be a general phenomenon in red giant and supergiant stars. 



5. CONCLUDING REMARKS 

In concluding, we summarize our present viewpoint: 

1. Water was discovered in several early M (super)giants nearly 40 years ago with 
Stratoscope II, but this discovery has been mis-interpreted and overlooked until 
recently. We regret that this important discovery did not provide major impact on 
the theory of cool stellar atmosphere during these 40 years. We certainly hope that 
such an unfortunate history should not be repeated with the ISO data. 

2. After 40 years, we confirmed H2O absorption bands in more than 30 red giant 
stars from K5III to M8III with the use of the ISO Data Archive. Thus ISO finally 
established that the presence of water is a general phenomenon in red giant stars 
including K and early M types. 

3. With ISO SWS, we detected water in emission in the M2 supergiant jj, Cep. This 
fact suggests that water should be in the outer atmosphere rather than in the pho- 
tosphere. Also, emission not only of H 2 but also of C0 2 and S0 2 are detected in 
late M giants with ISO and, by implication, water in all the red giant stars may also 
be originating in the outer atmosphere. 

4. We conclude that the presence of a rather warm molecular sphere - MOLsphere - 
may be a general feature in red giant and supergiant stars. Thus ISO revealed a 
novel picture of the atmosphere consisting not only of the photospher, chromosphere, 
and wind so far known but also of a new component - MOLsphere. Certainly, the 
ISO Data Archive should be an invaluable tool by which to explore the fundamental 
problem on the atmospheric structure of red giant and supergiant stars. 
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Figure 2: Observed high resolution SWS spectra (filled circles) are compared with the water 
spectra of Fig. 1 (solid lines). 




Figure 3: Observed low resolution SWS spectra (open circles) are compared with the water 
spectra of Fig. 1 (solid lines). 
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Figure 4: The observed N co \ values from the high and low resolution SWS spectra (filled 
and open circles, respectively) plotted against spectral types are compared with the predicted 
N co \ values based on the model photospheres (dashed line). 




Figure 5: Photospheric spectrum Fp (bottom) based on the spherically symmetric non-grey 
model photosphere (T e s = 3600 K, M* = 15M Q; R* = 650 R & ) is used as the boundary 
condition for solving radiative transfer in the MOLsphere. The resulting emergent spec- 
trum from the MOLsphere i*p + M (middle) is further diluted by the emission due to the 
optically thin dust envelope and Fp + m+d ( top ) is the final spectrum to be compared with 
the ISO spectrum f b s . The computations of Fp and -Fp+M o^e done with a resolution of 
R m 10 5 and the results are convolved with the slit function of SWS (R = 1600j. 
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Figure 6: A working hypothesis on the evolution of the atmospheric structure of red giants 
with Sp. Type. Stellar radius i?* is defined by tr oss ph 1 but photosphere extends to where 
tross ~ 0. Presently, self- consistent modelings are possible for the photosphere and interior 
( < t < oo), but not at all for other components in the outer atmosphere (—oo < t < 0). 



